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(February 8, 2008) 

We report on measurements by the E864 experiment at the BNL-AGS of the yields of light nuclei 
in collisions of 197 Au with beam momentum of 11.5 A GeV/c on targets of 20S Pb and 197 Pt. The 
yields are reported for nuclei with baryon number A=l up to A=7, and typically cover a rapidity 
range from y cm to y cm + 1 and a transverse momentum range of approximately 0.1 < Pt/A < 
0.5GeV/c. We calculate coalescence scale factors E>a from which we extract model dependent 
source dimensions and collective flow velocities. We also examine the dependences of the yields on 
baryon number, spin, and isospin of the produced nuclei. 
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I. INTRODUCTION 



Relativistic heavy ion collisions are believed to reach energy densities an order of magnitude greater than that of 
normal nuclear matter. These collisions allow the examination of the strong interaction in a novel environment as well 
as providing a possible doorway to new states of matter. In order to understand the dynamics of the collision system, 
one must use the only available tools-the species and momenta of the particles which exit the collision region. The use 
of emitted hadrons to probe the collision system is complicated because these hadrons rescatter many times as they 
traverse the collision region, and consequently lose some of their direct information about the earlier stages of the 
time evolution. However, the final space-time extent of the system at freeze-out (the time when strong interactions 
cease) and position-momentum correlations of the emitted particles contain much information about the entire time 
evolution. In principle, these carry information about the equation-of-state of the early collision region. 

In order to extract information about both the momentum and position distributions of the source at freeze-out, it is 
necessary to measure mult i- particle correlations. One widely used technique is Hanbury-Brown-Twiss interferometry 
for which the correlations between particles are due to quantum statistics. Another method is through the measured 
yields of light nuclei, which are formed by the coalescence of individual nucleons. 

Because of the violence of heavy ion collisions, it is highly improbable for a nuclear cluster near center-of-mass 
rapidity y cm =1.6 in a collision at these energies to be a fragment of the beam or target nucleus Q This would 
involve a cluster suffering a momentum loss of several GeV/c per nucleon that does not destroy the cluster, which 
is typically bound by only a few MeV per nucleon. These nuclei then arc formed by coalescence and so represent 
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correlations of several nucleons. As the mass of measured nuclei increases, of course, so does the number of particles 
involved in the correlation, and so does the sensitivity to features of the freeze-out distribution. 

In part due to the fragility of these states, the observed light nuclei are believed to be formed only near freeze-out 
of the collision system, at which time the mean free path of a bound cluster is long enough for it to escape without 
further collision. It is this notion that gives rise to a class of models of light nuclei production, the coalescence models 
(for example, @-|(|). In general, these models assume a phase space distribution of nucleons at freeze-out and impose 
some coalescence conditions on the freeze-out positions and momenta of the nucleons in order to calculate the yields 
of nuclei. These models differ both in their assumptions about the phase space profiles and coalescence conditions. 
Their differences are often characterized by their predictions of the invariant coalescence, or Ba, parameters which 
are defined as 

b {E ^ ] m 

where a nucleus with baryon number A and momentum P = pA is formed out of Z protons and N neutrons. 

In the simplest momentum space coalescence models, coalescence is assumed to take place between any nucleons with 
a small enough momentum difference. Early experimental results at the Bevalac with beam energies of » 500v4MeV 
and high energy proton induced reactions revealed Ba parameters which were approximately constant for these 
different collision systems ( contain useful summaries). The assumption of such simple models is that if the 
collision spatial volume is similar to size of the cluster, all nucleons whose momentum difference is less than a fixed 
value will fuse to form a nuclear cluster. Thus the experimentally observed constant Ba values may indicate collision 
volumes in these systems that are not substantially larger than the RMS size of a deuteron. 

In more advanced models, assuming a quantum mechanical sudden approximation j^] and using a density matrix 
formalism, accounting for both the positions and momenta of the nucleons 0] , the Ba parameters take on a relationship 
with the source volume V of Ba oc (l/V)^" 1 ^. Heavy-ion experimental results at higher energies at the AGS 
(sa lOAGeV) and CERN (w WOAGeV) revealed Ba values that decreased with beam energy. This observation was 
understood to be a sign of significant expansion in the collision volume before freeze-out. This larger source volume 
creates a situation where some nucleons with small relative momentum will have too large a spatial separation to 
coalesce, thus reducing Ba- 

The density matrix formalism [Q] assumes that although the collision volume can expand significantly, there is no 
correlation between the momentum and the position of a given particle. This assumption then leads to a prediction 
of no kinematic dependence of the Ba parameter. However, there is a great deal of evidence that collective motion is 
present leading to expansion of the collision volume and significant position-momentum correlations ]§| jj| . Although 
the overall expansion of the system tends to decrease Ba values by spatially isolating nucleons from each other, 
collective motion makes it more likely that nucleons that are spatially close together also have similar momenta, 
which to some extent works in the opposite direction by increasing Ba- Other coalescence models have made an effort 
to include the effects of both larger source volumes and collective motion, by including coalescence as a 'afterburner' 
in collision cascade models such as RQMD jlO],[| and by analytical calculations (TTJ fjlJ. 

Light nuclei production can also be calculated from thermal models |l3) 0,^1 which assume at least local thermal 
equilibrium of the system and thus that particle production (including, in this case, composite particles) is governed 
by a single temperature and chemical potentials. This gives rise to expressions for the Ba with the same (l/V)^ -1 ' 
dependence found in some coalescence models. Collective expansion can be included in these models ; it affects only 
the amount of energy available in local rest frames for particle production. 

There has been much experimental effort in the measurement of light nuclei in heavy ion reactions. Previous 
measurements at AGS energies, including clusters of A < 4, show values for Ba that are considerably lower than at 
Bevalac energies [lq| , indicating a much greater expansion of the system. This is also consistent with AGS results 
showing that the Ba become smaller with more central collisions and larger target nuclei |16| At CERN-SPS 



energies, production of secondary particles (chiefly pions) is several times as large as at the AGS, leading to a larger 
expansion before hadronic freeze-out. Coalescence of high mass clusters is thus much less probable, as indicated by 
yet lower values for the Ba as measured for example by experiment NA52 0. 

In this paper we describe and report the results of measurements by E864 of the yields of light nuclei in collisions 
of 197 Au with beam momentum of 11.5 A GeV/c on targets of 20S Pb and 197 Pt. The yields are measured for nuclei 
from baryon number A=l-7. In Section || we briefly describe the experimental apparatus and the analyses used to 



produce our final invariant multiplicities. In Section [II we report the results and compare them with measurements 



of other experiments where such measurements overlap. Finally in Section IV we examine trends in the data and 



discuss interpretations in the context of several different models of light nuclei production. 
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II. EXPERIMENT 864 



A. Apparatus 

Brookhaven AGS Experiment 864 is an open geometry, high data rate spectrometer which was chiefly designed 
to search for rarely produced objects in Au+Pb collisions. Figure |l| shows a schematic view of the experimental 
apparatus, a thorough description of which is given in Jl8| . Event centrality (impact parameter) is characterized 
by the charged particle multiplicity measured in an annular scintillator counter |19|| located approximately 10 cm 
downstream of the target, which subtends an angular region from 16.6° to 45° in azimuth when viewed from the 
target. The products of an interaction travel downstream through two dipole spectrometer magnets, Ml and M2. 
Charged particle identification is performed using information from the scintillator hodoscope walls (HI, H2, and H3) 
and the straw tube stations (S2 and S3). The hodoscopes walls each consist of 206 1cm thick scintillator slats placed 
vertically. They provide information about the charge, time-of-flight, and position of each charged particle hit, and 
this information is used to identify candidate charged particle tracks. These tracks are then rejected or confirmed and 
further refined by spatial hit information provided by the straw tube stations. Under the assumption that the track 
originates in the target, a rigidity is assigned to the track by a look-up table generated from a GEANT simulation 
of the experimental apparatus (using a technique described in [^o| as applied for the PHENIX experiment at RHIC). 
With information on rigidity, time-of-flight, and charge, a mass can then be assigned to the track, providing particle 
identification. A typical charge one mass distribution with a field of 0.45 Tesla in our spectrometer magnets is shown 
in Figure |2|. Mass resolutions of 3 to 7% RMS are typical for particles with velocity (3 < .985. 

At the downstream end of the apparatus is our hadronic calorimeter f2~H| ; an array of 754 towers, each measuring 
10cm x 10cm on the front face, each of which provides energy and time information. This is the essential piece of the 
apparatus in our analyses of yields of neutral particles, for which the tracking detectors serve only to provide a veto. 
In addition, the energy and timing information from each tower is used to provide the input for a level II high mass 
trigger (the Late-Energy Trigger or LET j22j), which provides an enhancement of approximately a factor of 50 in our 
searches for rare high mass states. 

Measurements reported in this paper are from a variety of experimental conditions, including different trigger 
conditions and different magnetic field settings in Ml and M2. The different data sets and their experimental 
conditions are listed in Table 1. Because of the large acceptance open geometry design of the experiment, the different 
data sets often have significant regions of overlap with one another, allowing a consistency check on the measurements- 



see for example the measurements of alpha particles introduced in Section III 



B. Data Analysis 

In order to measure the yield of a given species, a mass plot analogous to Figure || is made for each kinematic bin. 
The number of tracks which lie within its mass peak are determined. Background under each peak is then estimated, 
generally with fits to signal plus background, and subtracted away from this count. The invariant multiplicity in a 
given kinematic bin is then determined by correcting this number of raw counts for the geometric acceptance, trigger 
efficiency, charge cut efficiency, track quality cut efficiency, and detector efficiency. 

Geometric acceptances are generally 25% or lower as shown in Figure || for protons. Charge cut and track quality 
cut efficiencies are typically 90 to 95%. The three redundant charge measurements for each track allow easy calculation 
of the charge cut efficiency in each hodoscope simply by examining charge measurements made in each hodoscope 
against results from the other two; charge misidentification in each of the three hodoscopes is less than one percent, 
so fewer than one track in a million is assigned an incorrect charge. In general, the track quality cuts are determined 
by comparison with monte-carlo simulations. In cases where the efficiencies are particularly high, they are determined 
directly from the data. The total track detector efficiency is 85 to 90%; this is determined by excluding each detector 
in turn from the track-finding process. 

The trigger efficiency is significant only for those measurements which were made using the LET. This mass and 
momentum dependent efficiency ranges from approximately 40% to 90% for the measurements reported here. The 
LET efficiency is determined by two methods. The first, used mainly for slow, high mass states for which the efficiency 
is very high, is done simply using a monte-carlo simulation of the shower generated by the object and knowledge of 
the LET look-up table in each tower. The second method is from the equation clet — Nlet / (Nlet + R X N non LET) 
where R is the rejection factor provided by the LET (i.e. number of events which fire the LET divided by the total 
number of events). Nlet and N non LET are the numbers of particles of interest which do and do not produce LET 
triggers in LET triggered events, respectively. 
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Sources of systematic error that we have quantified include possible error in the determination of the efficiencies 
listed above as well as error in background subtraction (particularly relevant for the deuterons and alpha particles 
at their highest rapidities). Also examined were the effects of changing the assumed input distribution for each 
particle species in the determination of geometrical acceptances and efficiencies and effects of possible differences in 
the magnetic field with the field maps that were used in reconstruction of tracks. 

Overall errors are generally dominated by systematics, particularly for the lighter states. Statistical errors can be 
significant for the heavier states, particularly in the determination of LET efficiencies in which the number of particles 
of interest which do not fire the LET generally has the largest statistical error. 



III. RESULTS 



Measurements of invariant multiplicities for protons, neutrons, deuterons, s He, and 4 He are shown in Figures Q 
through ||, and the values of the data shown in these figures are listed in tables in the Appendix. Figure ^ displays 
proton invariant multiplicities for three different bins of collision centrality. Figure shows only proton yields from 
10% most central Au+Pb collisions along with neutron multiplicities measured by E864 (see Reference p3[) for 
comparison. Figures || and display deuteron and 3 He invariant multiplicities for the same centrality bins used for 
the proton measurements. Tritons, 4 He, 6 He, 6 Li, 7 Li, and 7 Be are measured by E864 only in 10% most central 
collisions; yields for tritons and alpha particles are shown in Figure || and Figure g, respectively, while yields for the 
heavier nuclei are listed in tables in the Appendix. Added detail concerning most of these measurements may be 
found in the Ph.D. theses listed in Reference plfl . 



A. Contributions from Hyperon Decays 



For comparison to other experimental results and calculations of coalescence parameters, it is important to quantify 
the contribution to proton yields that is made by protons which come from decays of hyperon states, which to E864 
are indistinguishable from primordial protons. There are three dominant hyperon decays which produce protons: 
A — > p + n~ , S° ^A + 7^p + 7r~+7, and E + — > p + ir°. The contributions from these decays were evaluated 
using a GEANT simulation of the experiment with an input distribution taken from measurements of E891 for the 
Lambda |^5| and an input distribution from RQMDv2.3 |l(]] for the Sigma. From this simulation it was determined 
that protons from hyperon decays account for approximately 12% of the measured yields of protons with only a slight 
kinematic dependence. The proton and neutron yields from E864 shown in Figures ^ and |^ have not been corrected 
for hyperon feed-down (nor have the values listed in the appendix tables). 



B. Comparisons with Other Experimental Results 



Figure |ic| shows a comparison of the light nuclei measurements from AGS experiments E864, E877 pq] , and E878 
p6[ . Because of the different beam momenta of the experiments (10.8 A GeV/c in E878), the yields are plotted versus 
beam normalized rapidity. The yields shown for E864 and E877 are average yields at approximately pt/A =150 
MeV/c, while the E878 yields are measurements at px — 0. Other caveats to the comparisons of the yields shown in 
Figure 10 are noted in the figure caption. 

Proton yields measured by E864 are clearly higher than measurements by E878. Some of this difference can be 
attributed to protons which are feed-down from hyperon decay. The acceptance of E878 for these feed-down protons 
is only about 10% of what it is for primordial prot ons, w hile in E864 the two acceptances are nearly the same. When 
this difference is taken into account (see Section III A) and the E864 yields of primordial protons are lowered by 
approximately 12%, the results of the two experiments are different by approximately 25% at midrapidity; this lies 
within the range of systematic errors for the two experiments. At higher rapidities, the agreement is better. 

Comparison of the three experiments' measurements of deuterons and tritons shows close agreement and the mea- 
surements of E864 and E878 of alpha particle yields also agree within errors. 



IV. DISCUSSION 



A. General Trends of the Spectra 
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1. Transverse Dependence of Yields 



E864 has sufficient coverage in transverse momentum for us to extract measurements of inverse slope parameters of 
the different light nuclei yields in the rapidity range 2.2 < y < 2 .4. Show n in Figure |ll| are yields of protons, deuterons, 
3 He, and He as a function of transverse mass tut — mo = \/p\ + itiq — too in this rapidity slice. Overlayed on the 
measurements are fits of each species to a Boltzmann distribution in transverse mass, from each of which we extract 
an inverse slope parameter, T, as noted in Figure O. The fits from which we extract T are linear fits to the log of 
the invariant multiplicities divided by tot; this is the same fitting method used for determining the slope parameters 
for neutrons in [|23|| . The x 2 values are less than one per degree of freedom for all these fits. 

These slope parameters, as well as those for neutrons in this same rapidity range, are displayed in Figure ^ as a 
function of mass number. Polleri et. al. j[2| have demonstrated the sensitivity of these trends in the inverse slopes 
to the density and velocity profiles of the nucleons at the time when coalescence occurs. To make this point, they 
have performed calculations of the behaviour of these trends for different assumptions about the source distributions. 
Two of these assumptions give rise to the two curves shown overlayed on the data in Figure |l2|. The first has a 'box' 
spatial profile and a linear velocity profile, and the second has a Gaussian spatial distribution and a velocity profile 
v(tt) oc r^ 2 \ In determining the curves shown in Figure 12 we have fit the data to these two different functional 
forms with the numerical constraint of T=100 MeV for zero mass. Neither of these sets of model parameters provides 
an adequate description of the data. These curves are meant only to illustrate the sensitivity of these measurements: 
clearly, the sensitivity to the differences in these assumptions increases with increasing mass of the measured nuclei. 

Shown in Figure |l3| are these same trends in light nuclei slope parameters including centralities other than the 10% 
most central collisions. For the more peripheral events, the trends are consistent with a linear dependence of slope 
parameter on mass (matching the calculation in Reference fl2|| including box density profiles). Only for the most 
central events is there a clear rollover in the slope as a function of mass. 



2. Longitudinal Dependence of Yields 

In order to examine the rapidity dependence of the yields of light nuclei, we observe the trends in multiplicities in the 
Pt/A range from 100 to 200 MeV/c. (Our transverse coverage is not sufficient to integrate in transverse momentum 
for a full measurement of dN/dy over this entire rapidity range.) In Figure |lj we show the invariant multiplicities of 
deuterons for this low px range for three different centralities as a function of rapidity. The yields are concave as a 
function of rapidity (i.e. they are lowest at center-of-mass and increase toward beam or target rapidity) and become 
more concave for the more peripheral events. We can parameterize this concavity by fitting each set of yields to a 



quadratic a + b(y — y cm ) 2 . These fits are shown overlayed on the data in Figure 14. The ratio of coefficients, b/a, from 
this parameterization serves as a measure of the relative concavity of a species' spectrum as a function of rapidity at 
lowp T |27j]. 

Values of this ratio b/a are plotted in Figure [l5| as a function of mass number A for protons, deuterons and 3 He in 
three different collision centralities. We observe that the concavity of a spectrum increases with the mass of a species 
and with collision centrality. This can be understood as some form of collective motion in the longitudinal direction, 
either expansion or incomplete stopping, which increases in peripheral events. It also is consistent with a rapidity 
dependent transverse expansion which pushes the nuclei out of the transverse momentum range measured here. 



B. Ba parameters 



The trends noted in Section IV A will also be observable in a study of the behaviour of the Ba parameters. For 
simplicity, we can from our neutron measurements [ p3[ characterize the neutron spectrum as a factor of 1.19±.08 
greater than the proton spectrum (feed-down from hyperons is taken into account in determining this ratio), and so 
we evaluate Ba as 

( F d 3 N A \ 

= Z7ZJpn=77- ( 2 ) 



(1.19) N (E 



dp 3 



for a nucleus of baryon number A with N neutrons and Z protons. Again, all invariant yields are evaluated at a 
common velocity. 
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In Figures |16| and |l7| we show measurements of B2 and B3 as a function of rapidity and transverse momentum. 
Because nucleons which are feed-down from weak decays are not available as nucleons for coalescence, we have 
subtracted the contribution to proton and neutron invariant yields from hyperon feed down for the calculation of Ba- 

The values of the Ba in Figures ^ and ^ are clearly not momentum independent as is assumed in many early 
coalescence models ( 0] and ||, for example). This is expected given the large amount of evidence for collective 
flow in heavy- ion collisions [^) . The values of the Ba parameters seem to increase slightly with increasing transverse 
momentum and there is a clear increase away from center-of-mass rapidity. Both of these increases in Ba away 
from center-of-mass momentum are consistent with expectations from an expanding source, although the longitudinal 
dependence can also be interpreted as a sign of incomplete stopping. Indeed the fact that the invariant yields of 
antiprotons near pt=0 are strongly peaked near center-of-mass rapidity as measured by E864 |^8| while the proton 
yields are essentially flat may be taken as further evidence of incomplete stopping. 



C. Source Size Calculations 



The Ba parameters can be related through coalescence models to source sizes. Following the model of Sato and 
Yazaki Q which uses a density matrix representation of the source distribution and projects it onto a representation 
of the dcuteron wave function, we can relate the Ba parameters to the RMS source radius through 

B A = ( 2 ^)^^^)^^ (3) 
2 A m\ A L > va + v 

with va the size parameter for a cluster with baryon number A and spin J a (m represents the nucleon mass). This 
model assumes the absence of collective motion of the nucleons, and therefore a radius (given by Rrms = \/3/2v) 
which is independent of momentum. We evaluate the source size using B2, -B3, B4, Bq, and B-j values from our 
measurements nearest the collision center-of-mass (again assuming a neutron to proton ratio of 1.19) and list the 
results in Table ||. We have used values for the va from Q for A — 2,3 and 4 and following (29) have done a 
polynomial extrapolation to determine vq and V7. This extrapolation for va may be suspect particularly for the halo 
nucleus 6 He, but the final value for Rrms is quite insensitive to the value for vq. The extracted radius for deuterons 
A=2 is considerably larger than the initial size of the colliding nuclei, and the radius parameters decrease for clusters 
of increasing mass. 

A model by Scheibl and Heinz JTT[ which includes the effect of collective flow in a density matrix prescription for 
coalescence leads to an expression for source dimensions: 

37T 3 / 2 < C d > 2(mr-m)(^-^) 
2 2m T Rl(m T )R ll (m T ) e * U 

where Rj_ and Um are the transverse and longitudinal dimensions of the fraction of the source which contributes to 
deuteron emission (comparable to the radius parameters extracted in the YKP parameterization of HBT interferom- 
etry) , < Cd > is a quantum mechanical correction factor for the finite size of the deuteron which is evaluated by the 
authors under various assumptions about the source, and T* and TjJ are the inverse slope parameters for protons and 
deuterons. Equation^ as written assumes a box density profile for the source; a Gaussian profile would result in the 
absence of the final exponential factor. Plugging in our results for B2, we can extract values for (i?5_(wr)-R|| (^t)) 1 ^ 3 



which are shown in Figure 18 as the solid circles. For the calculations shown here we have used 0.75 for < Cd > and 
the values for and T* as measured at rapidity 2.3. 

Shown also in Figure [l^ as hollow circles are the results of source size calculations with the similar fragment 
coalescence model of Llope et. al. |^9| via the equation 

3 3/2 (2J e + l) m c — 

tt r - = ' — — ; T- : tttITt 



(2J a + l)(2J b + 1) m a m b Ec ^ 



a-jVq 771 < 



which relates the effective source radius R in the frame of a cluster c which may be formed through the coalescence 
of smaller clusters a and b. (Note that the radius parameters R shown in Figure [l8] are meant to describe sources of 
the form p(r) oc exp(—r 2 /2R 2 ) and so correspond to RMS radii of Rrms = v3R.) 

We observe in Figure [fq that there is a decrease in source size with increasing distance away from the center-of-mass 
(again, as expected for a source with radial expansion), but our measurements do not give a clear picture concerning 
scaling of the source size with transverse mass ]2^ ] such as has been noted in results for sizes from HBT two-particle 
correlations and in measurements of B2 at the CERN SPS |30|. 
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D. Comparison with RQMD 



We can also compare our results with predictions from the cascade model RQMD version 2.3. The complex 
many-body processes by which light nuclei are formed are not included in RQMD, rather an afterburner 31 1 is used 
to calculate the coalescence of these states based upon an input of the positions and momenta of the nucleons at 
freeze-out. This model then explicitly includes the position-momentum correlations due to expansion, etc. that are 
present in RQMD. 

Figure [l9| displays the E864 measurements along with predictions from RQMD with the coalescence afterburner 
for comparison over the transverse momentum range 0.1 < pr/A < 0.2 GeV/c as a function of rapidity. For 
comparison, RQMDv2.3 was run under two different conditions, one including the effect of repulsive mean-field 
potentials (potentials mode) and one not (cascade mode). We note that there is an increasing disagreement with 
increasing mass as noted previously in Reference |l6|. Calculations of 4 He production in cascade mode are at a level 
of approximately 100 lower than our measurements; with potentials on, the discrepancy is still larger by about a factor 
of two. The level of disagreement varies somewhat in with transverse momentum (the slopes are in fact generally 
better predicted with potentials on than off) but particularly for the heavier states the change is slight compared to 
the overall level of disagreement. 



E. Scaling of Yields Versus Mass 

Shown in Figure |2^ are the invariant yields in a small kinematic region at or near y — 1.9, pr < 300MeV. Over 
ten orders of magnitude, the yields in this kinematic bin fit very closely to an exponential dependence with a penalty 
factor of approximately 48 for each nucleon added (see |32| and references therein for a discussion of such exponential 
behavior of cluster yields at lower collision energies). Of course, as we have seen above, the yields for different species 
have different kinematic dependences due to collective motion, and so this should not be taken as a penalty factor 
which governs the integrated yields of the various species, which from the kinematic dependences of the Ba parameters 



discussed in Section IV B can clearly be different. 

We can make a crude estimate of the difference in these two penalty factors by parameterizing the rapidity and 
mass dependences of the inverse slope parameters as follows: from our neutron measurements |23[] which extend up 
to beam rapidity, we can roughly parameterize the distribution of T in rapidity as a Gaussian with a width of 1.1 
units. From Figure we also can make a rough parameterization of the mass dependence of the inverse slope as 
T oc (A + 1.0) With these two parameterizations and our penalty factor of 48 at y = 1.9, pt < 300MeV, we calculate 
a penalty factor of 25 in the integrated yields for the addition of an extra baryon to a coalesced state. This should be 
considered as a lower limit on such a penalty factor since as previously noted the inverse slope par ameter rolls over 
as a function of mass rather than following the linear dependence we have used in this estimate [B3fl. 

We can also estimate this penalty factor in overall yields by following section VI of Reference pll|. Equation 6.10 
in this reference allows us effectively to make an estimate of the difference between the penalty factor near px = 
and the overall penalty factor in two limiting cases: a static, homogeneous fireball (which would translate our penalty 
factor of 48 at pr =0 to an overall penalty factor of 72) and a rapidly expanding system (which would result in an 
overall penalty factor of 39). 



F. Spin and Isospin Dependences of Yields 

In Figure |l|, we display three ratios as a function of transverse momentum in the rapidity ranges y=l. 8-2.0 and 
y=2.0-2.2. The three ratios are the ratios of invariant yield of neutrons over invariant yield of protons, the ratio of 
yield of tritons over yield of 3 He and the ratio of 6 He to 6 Li. The n/p and t/ 3 He ratios are consistent with a value 
of approximately 1.2 (in Reference p3| , we extract values of 1.19 ± .08 and 1.23 ± .04, respectively, for the two ratios 
in the range 1.6 < y < 2.4.). In contrast, the 6 He/ 6 Li is much nearer to a value of 0.3. p, n, 3 He and t are all spin 
J=l/2 states while 6 He is spin and 6 Li is spin 1. We take this as evidence that the yields scale as the degeneracy 
factor 2J+1 which is commonly predicted in thermal and coalescence models. 

With the dependences upon mass number, isospin and spin divided away one can examine the yields for other 
dependences. This topic, including a possible dependence on binding energy per nucleon with an inverse slope 
parameter dependence of a few MeV, is discussed in Reference p4j . 
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V. SUMMARY 



We have shown results of measurements of light nuclei from A=\ to A=7. The increase with mass of light nuclei 
inverse slope parameters appears to roll over at approximately A— 3 in central events but not necessarily in more 
peripheral events. Also we have shown that the yields near px—0 are concave as a function of rapidity and that this 
relative concavity increases in more peripheral events and in higher mass nuclei, consistent with both radial expansion 
and incomplete stopping. These trends are also evident from our observations of the kinematic dependences of the 
Ba parameters. From these parameters we have extracted source dimensions from various models. Efforts to extract 
more quantitative information about the source from these measurements using the cascade model RQMD with a 
coalescence afterburner were unsuccessful as predictions of the model differ from our results by an amount that 
increases with mass and reaches a level of 100 or more by A=A. 

We have also examined the overall scaling of the yields up to A— 7, extracting a penalty factor of about 48 to add 
a nucleon to a coalesced state near midrapidity at low transverse momentum. This likely translates into a somewhat 
smaller penalty factor in overall yields for the addition of a nucleon, but we have argued that this is unlikely to differ 
by as much as a factor of two from our measured penalty of 48. 
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Shown in Table III through Table XVI arc results of measurements by E864 of invariant multiplicities of light nuclei 
from A=l through A=7 in 10% most central Au+Pb collisions. Results for less central data are also listed for protons, 
dcutcrons and 3 He. 
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FIG. 1. The E864 spectrometer in plan and elevation views, showing the dipole magnets (Ml and M2), hodoscopes (HI, H2, 
and H3), straw tube arrays (S2 and S3) and hadronic calorimeter (CAL). The vacuum chamber is not shown in the plan view. 
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FIG. 4. Panel (a) displays invariant yields for protons in 10% most central Au+Pb collisions. Panels (b) and (c) display the 
same for 10-38% and 38-66% central collisions, respectively. 
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extracted radii from the cluster coalescence model of Llope et. al. 
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FIG. 19. Yields measured by E864 near pr = (in units of <? /GeV 2 ) as a function of rapidity for protons (circles), deuterons 
(squares), 3 He nuclei (stars) and alpha particles (triangles). For comparison, the predictions of the model RQMDv2.3 with 
a coalescence afterburner are shown for the four species both for RQMD in cascade mode (solid histogram lines) and with 
mean-field potentials (dashed histogram lines). Hollow symbols represent reflections of E864 data points about midrapidity. 
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FIG. 20. Mass dependence of invariant yields of light nuclei from A=l up to A=7. Yields are measured at y=1.9, pr/A < 
300 MeV/c except where otherwise noted in the plot. 
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FIG. 21. n/p, t/ 3 He, and 6 He/ 6 Li ratios as a function of transverse momentum in two different rapidity slices. 
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TABLE I. Summary of experimental conditions under which the light nuclei measurements were made. Events refers to the 
number (in millions) of sampled events in a given data set. 



Year 


Field (T) 


Trigger 


Events (M) 


Species 


Rapidity 


Centrality 




1994 


+.75 


MULT 


24 


p,d, 3 He*He 


1-2 < y < 2.0 


0-10% 




1995 


+.45 


MULT 


6 


p,d, 6 He 


1-2 < y < 2.4 


0-10,10-38,38-66% 




1996 


+1.5 


MULT 


7 


n 


1-6 < y < 3.2 


0-10.10-38,38-66% 




1996 


+1.5 


MULT+LET 


13000 


"He, b He, a Li, ' Li,' Be 


1-6 < y < 2.2 


10% 




1998 


+.45 


MULT+LET 


2000 


4 He 


1-4 < y < 2.4 


10% 





TABLE II. Ba parameters in units of (GeV 2 /c 3 )^ A 1 ' near collision center-of-mass and source sizes extracted from them 
from the model of Sato and Yazaki. 



A 


species 


B A 


RnMs(fm) 


y,pr 


2 


d 


8.7+1.4 x 10~ ; ' 


14.8+ .7 


1-4 < y < 1.8, £f < 200MeV 


3 


*He 


5.1 ± 1.1 x 10"' 


12.2+ .5 


1-4 < y < 1.8, 2f < 2QQMeV 


4 


"He 


4.9 ± 1.3 x 10~ 1U 


10.1+ .3 


1-4 < y < 1.8, £f < 200MeV 


6 


b He 


1.9 + 0.8 x 10~ lb 


9.9+ .3 


1-6 < y < 1.8, £f < 170MeV 


7 


'Be 


1.4 + 0.7 x 10" is 


9.2+ .3 


1-6 < y < 2.2, £f < 250MeV 



TABLE III. Invariant yields for protons in 10% most central Au+Pb collisions in units of c 2 /GeV 2 Errors listed are systematic 
and statistical combined in quadrature. 



Rapidity 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


2.4-2.6 


jj T (MeV/c) 














25-50 


17.9+ 3.5 


22.4+ 9.6 










50-75 


28.0+ 3.5 


25.0+ 2.2 


25.2+ 3.5 


14.9+ 5.1 






75-100 


26.0+ 2.2 


27.5+ 2.4 


25.8+ 2.2 


26.2+ 2.5 


32.1+ 6.1 




100-125 


25.2+ 3.7 


25.9+ 1.9 


27.1+ 2.0 


27.0+ 2.5 


23.3+ 2.4 


20.8+ 4.3 


125-150 




26.2+ 3.0 


25.4+ 1.7 


26.2+ 1.9 


24.2+ 1.9 


25.2+ 2.3 


150-175 




23.8+ 6.4 


26.4+ 2.0 


27.0+ 1.9 


26.4+ 1.9 


26.1+ 2.1 


175-200 






24.6+ 1.9 


24.3+ 1.8 


26.3+ 1.9 


26.5+ 2.1 


200-225 






23.4+ 4.1 


25.4+ 1.9 


25.4+ 1.8 


26.0+ 1.9 


225-250 






21.1+ 2.2 


24.1+ 1.9 


24.9+ 1.7 


26.2+ 1.8 


250-275 






23.8+ 4.6 


23.1+ 1.7 


23.2+ 1.7 


25.5+ 1.9 


275-300 








21.9+ 1.6 


23.5+ 1.9 


24.5+ 1.7 


300-325 








21.1+ 2.0 


22.3+ 1.6 


22.0+ 1.6 


325-350 








20.0+ 1.6 


22.6+ 1.6 


20.7+ 1.6 


350-375 








17.6+ 2.7 


22.8+ 1.7 


19.7+ 1.5 


375-400 








15.1+ 2.0 


21.3+ 1.6 


19.7+ 1.5 


400-425 










21.6+ 1.6 


17.8+ 1.5 


425-450 










17.8+ 1.3 


16.3+ 1.3 


450-475 










17.8+ 1.5 


16.7+ 1.4 


475-500 










17.0+ 2.5 


16.3+ 1.4 


500-525 










15.9+ 1.4 


15.3+ 1.2 


525-550 












14.4+ 1.2 


550-575 












13.5+ 1.1 


575-600 












13.5+ 1.1 


600-625 












12.0+ 1.0 


625-650 












10.0+ 0.8 


650-675 












9.2+ 0.9 
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TABLE IV. Invariant yields for protons in 10 to 38% most central Au+Pb collisions in units of c 2 /GeV 2 Errors listed are 
systematic and statistical combined in quadrature. 



Rapidity 


1 A 1 R 
1.4-1.0 


l.O-l.o 


i q o n 
l.o-z.U 


o n o o 
z.u-z.z 


O O O A 
Z.Z-Z.4 


p T (MeV/c) 












25-50 


14.3± 2.9 










50-75 


17.9± 2.3 


17.7± 1.6 


16.6± 2.3 


19. 1± 7.4 




75-100 


16. 5± 1.4 


16. 2± 1.4 


16.7± 1.5 


17.7± 1.9 


20.0± 4.2 


100-125 


16.8± 2.5 


16.8± 1.3 


17.3± 1.3 


19. 2± 1.8 


17.4± 1.8 


125-150 




16.6± 1.9 


16.8± 1.1 


16.6± 1.2 


17.2± 1.4 


150-175 




14.7± 4.0 


17.2± 1.3 


17.1± 1.2 


18.6± 1.4 


175-200 






14.4± 1.1 


16.3± 1.2 


18.2± 1.4 


200-225 






14.9± 2.6 


16. 5± 1.2 


17.1± 1.2 


225-250 






13.8± 1.5 


15. 1± 1.2 


17.4± 1.2 


250-275 






13.2± 2.6 


15.3± 1.2 


15. 5± 1.1 


275-300 








14.5± 1.1 


14.5± 1.2 


300-325 








13.0± 1.2 


14.5± 1.1 


325-350 








12. 1± 1.0 


14.9± 1.1 


350-375 








11. 3± 1.8 


14. 5± 1.1 


375-400 








10. 2± 1.4 


14.2± 1.1 


400-425 










13.8± 1.0 


425-450 










12.1± 0.9 


450-475 










10.4± 0.9 


475-500 










11. 1± 1.6 


500-525 










10.3± 0.9 



TABLE V. Invariant yields for protons in 38 to 66% most central Au+Pb collisions in units of c 2 /GeV 2 Errors listed are 
systematic and statistical combined in quadrature. 



Rapidity 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


p T (MeV/c) 












25-50 


4.4± 0.9 










50-75 


6.9± 0.9 


7.1± 0.7 


7.5± 1.1 


7.5± 2.9 




75-100 


6.9± 0.7 


7.2± 0.7 


7.4± 0.7 


7.6± 0.8 


9.2± 2.1 


100-125 


6.6± 1.0 


7.0± 0.6 


7.3± 0.6 


8.1± 0.8 


8.0± 0.9 


125-150 




7.3± 0.8 


7.0± 0.5 


7.4± 0.6 


8.3± 0.7 


150-175 




5.6± 1.5 


7.1± 0.6 


7.6± 0.6 


8.9± 0.7 


175-200 






5.7± 0.5 


6.8± 0.5 


8.5± 0.7 


200-225 






6.3± 1.1 


7.1± 0.5 


8.1± 0.6 


225-250 






5.1± 0.6 


6.6± 0.6 


8.0± 0.6 


250-275 






5.8± 1.2 


6.5± 0.5 


7.1± 0.5 


275-300 








6.6± 0.5 


7.0± 0.6 


300-325 








5.7± 0.6 


6.8± 0.5 


325-350 








5.3± 0.5 


7.1± 0.5 


350-375 








4.7± 0.7 


6.5± 0.5 


375-400 








4.3± 0.6 


6.3± 0.5 


400-425 










6.0± 0.5 


425-450 










5.2± 0.4 


450-475 










4.6± 0.4 


475-500 










4.8± 0.7 


500-525 










3.9± 0.4 
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TABLE VI. Invariant yields for deuterons in 10% most central Au+Pb collisions in units of 10~ 2 c 2 /GeV 2 in bins with low 
transverse momentum (see also the following table for the high transverse momentum data). Errors listed are systematic and 
statistical combined in quadrature. 



Rctpidity 


i n i 9 

1 . U- 1 . z 


19 1/1 
_L .Z- 1 


1 A 1 (\ 


i ^ i a 

l.O-l.o 


1 OOfi 


9 9 9 
z.u-z.z 


9 9 9 /1 
Z.Z-Z. 1 * 


9/1 

Z . Z . u 


p T (MeV/c) 


















37.5 


74. 7± 35.7 
















62.5 


76. 4± 8.8 


53. 7± 11.0 














87.5 


71. 2± 6.4 


63. 1± 5.8 


48. 2± 6.4 


75. 9± 26.7 










112.5 


69. 3± 6.2 


61. 7± 5.3 


52. 6± 5.8 


56. 1± 9.3 


43. 0± 9.7 








137.5 


72. 5± 8.6 


66. 0± 5.9 


53. 9± 5.9 


56. 6± 5.2 


54. 4± 5.6 


68. 1± 14.0 






162.5 


86. 6± 16.6 


59. 3± 4.7 


53. 2± 4.9 


54. 4± 5.4 


56. 4± 5.2 


70. 4± 6.9 


59. 1± 18.2 




187.5 




65. 1± 5.3 


58. 0± 4.6 


56. 5± 5.0 


53. 5± 4.4 


64. 5± 6.0 


70. 6± 8.6 




212.5 




67.5± 7.2 


50. 1± 4.3 


54.5± 5.2 


56.6± 4.9 


60.5± 4.7 


71. 0± 7.0 


89. 7± 13.7 


237.5 




56. 7± 8.4 


57. 3± 5.3 


59. 4± 4.6 


58. 0± 4.5 


61. 6± 4.7 


66. 5± 6.2 


91. 0± 9.6 


262.5 




67. 3± 14.3 


57. 2± 5.0 


59. 4± 4.7 


58. 6± 4.2 


65. 1± 4.9 


72. 9± 5.9 


98. 7± 8.2 


287.5 






56.3± 5.0 


52. 1± 5.2 


58.0± 4.2 


61. 1± 4.3 


74.0± 5.5 


86.3± 6.8 


312.5 






50. 8± 4.7 


51. 0± 4.5 


54. 5± 4.6 


64. 3± 4.3 


75. 2± 5.5 


95. 1± 8.4 


337.5 






46. 3± 5.7 


55. 1± 4.5 


56. 1± 4.8 


64. 6± 4.4 


72. 2± 5.1 


89. 8± 9.3 








co 74- a 7 




KA 84- c: c 


«fi f?4 A A 




ok a J- q 1 


387.5 






44.1± 11.3 


53.7± 4.5 


54.2± 4.5 


68.2± 4.9 


70.3± 5.1 


87.2± 8.8 


412.5 








47.8± 4.3 


55.1± 5.0 


61.8± 4.7 


71.3± 5.0 


82.0± 7.8 


437.5 








51.3± 4.4 


56.2± 4.6 


62.0± 5.4 


75.9± 5.2 


88.1± 9.1 


462.5 








48.2± 5.0 


56.5± 4.3 


58.3± 5.1 


66.7± 4.8 


85.7± 8.6 


487.5 








50.8± 7.9 


57.2± 4.4 


57.6± 4.5 


68.2± 5.2 


82.0± 8.8 


512.5 








46.1± 6.4 


55. 0± 4.3 


64.4± 5.3 


65. 2± 5.2 


80.3± 7.1 


537.5 










54.1± 4.3 


61. 1± 5.7 


59.5± 4.7 


75.0± 8.0 
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TABLE VII. Invariant yields for deuterons in 10% most central Au+Pb collisions in units of 10~ 2 c 2 /GeV 2 in bins with high 
transverse momentum (see also the previous table for the low transverse momentum data). Errors listed are systematic and 
statistical combined in quadrature. 



Rapidity 


i n i o 
1 .u- i .z 


19 14 
1 .Z- 1 .4 


1 A 1 (\ 
1.4- 1 .U 


1 (\ 1 8 


1 o 9(1 
l.O-Z.U 


9 9 9 
z.u-z.z 


9 9 9 /1 
Z.Z-Z.4 


9 A 9 f\ 
Z.4-Z.U 


p T (MeV/c) 


















512.5 








46. 1± 6.4 


55. 0± 4.3 


64. 4± 5.3 


65. 2± 5.2 


80. 3± 7.1 


537.5 










54. 1± 4.3 


61. 1± 5.7 


59. 5± 4.7 


75. 0± 8.0 


562.5 










51. 7± 4.3 


61. 2± 7.7 


64. 8± 5.6 


77.1± 8.6 


587.5 










50. 2± 4.4 


64. 4± 4.9 


62. 1± 5.3 


71. 1± 7.4 


612.5 










57.6± 4.9 


63. 5± 6.1 


60. 4± 5.2 


65. 9± 7.3 


637.5 










56. 3± 8.1 


55. 9± 4.2 


63. 6± 5.6 


66. 6± 8.4 


662.5 










43. 6± 4.4 


61. 8± 5.6 


58. 8± 5.2 


63. 7± 7.9 


687.5 










43.0± 10.0 


55. 7± 4.4 


61. 8± 5.8 


65. 7± 9.6 


712.5 












59. 5± 4.9 


58. 1± 5.1 


58. 4± 8.2 


737.5 












57. 4± 4.6 


55. 8± 4.6 


58. 0± 8.2 


762.5 












60. 2± 4.5 


53.1± 4.5 


55.4± 8.4 


787.5 












52. 6± 4.0 


53. 8± 4.2 


58. 6± 9.2 


812.5 












50. 7± 6.9 


54. 3± 4.3 


55. 8± 8.2 


837.5 












49. 3± 4.2 


55. 1± 4.7 


47.6± 6.9 


862.5 












43. 3± 4.9 


52. 8± 4.4 


45. 5± 5.8 


887.5 












51. 3± 7.0 


48. 4± 4.1 


50. 4± 6.7 


912.5 












47. 0± 8.4 


50. 5± 4.4 


44. 4± 6.3 


937.5 














49. 4± 4.3 


42. 6± 5.3 


962.5 














45. 7± 3.8 


46. 2± 6.9 


987.5 














44. 1± 3.9 


47.6± 7.5 


1U1Z.0 














AQ 94- A O 
4o.Z± 4.U 


A A 04- R O 

44. Z± o.y 


1037.5 














42. 2± 4.0 


40.7± 6.6 


1062.5 














38.1± 3.6 


40.6± 6.5 


1087.5 














37.0± 4.6 


40.9± 5.3 


1112.5 














39.9± 4.1 


40.7± 5.9 


1137.5 














37.6± 4.3 


36.8± 4.8 


1162.5 
















32.0± 4.7 


1187.5 
















31.9± 5.1 
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TABLE VIII. Invariant yields for deuterons in 10% to 38%most central Au+Pb collisions in units of 10 2 c 2 / GeV 2 . Errors 
listed are systematic and statistical combined in quadrature. 



Rapidity 


l.U-l.Z 


10 1/1 
l.z-1.4 


1 A 1 R 
1.4-1.0 


1 R 1 S 
1.0-1.O 


icon 
l.o-z.U 


o n o o 
z.u-z.z 


Z.Z-Z.4 


p T {MeV/c) 
















25 
















75 


48. 2± 5.2 


42. 5± 5.5 












125 


52.4± 5.3 


40.6± 3.8 


36.0± 4.2 


36.1± 4.7 


35. 2± 4.9 






175 




44.9± 4.5 


33.6± 3.6 


36.9± 3.5 


38. 1± 3.5 


42. 7± 4.6 




225 




42. 8± 5.4 


34.6± 3.3 


35. 2± 3.0 


36.4± 3.2 


38. 5± 4.0 


52. 3± 4.9 


275 






38. 2± 3.7 


34.4± 3.3 


38. 5± 2.9 


45.3± 3.5 


57.7± 4.3 


325 






33.4± 3.7 


37.7± 3.2 


35. 5± 3.6 


41.3± 3.4 


52. 1± 4.7 


375 






31.4± 6.9 


33.0± 3.1 


36. 2± 3.4 


42.3± 3.2 


56.4± 4.1 


425 








29.9± 2.9 


37.8± 3.7 


40.9± 3.4 


53.0± 3.9 


475 








29.4± 4.1 


35.9± 3.1 


37.0± 3.3 


46.8± 3.9 


525 








25. 0± 5.1 


35. 7± 3.0 


36.9± 4.6 


45. 2± 4.1 


575 










31. 5± 3.1 


41. 2± 3.5 


43.3± 4.2 


625 










33.6± 3.9 


34.2± 4.0 


43.3± 3.8 


675 










25.8± 4.3 


39.3± 3.3 


38.7± 3.5 


725 












33.6± 3.2 


42. 2± 3.7 


775 












33. 1± 3.6 


36.3± 3.1 


825 












27. 1± 3.5 


36.1± 3.9 


8(0 












oo o_i_ o a 
ze.3± 3.6 


Q A 0^_ O C 

34. l± 3.0 


925 












29.6± 4.2 


30.3± 3.9 


975 














25. 7± 2.4 


1025 














27.6± 2.5 


1075 














24.0± 2.5 


1125 














22.8± 3.1 


1175 














17.6± 2.8 
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TABLE IX. Invariant yields for deuterons in 38% to 66%most central Au+Pb collisions in units of 10 2 c 2 / GeV 2 . Errors 
listed are systematic and statistical combined in quadrature. 



Rapidity 


l.U-l.Z 


10 1/1 
l.z-1.4 


1 A 1 R 
1.4-1.0 


1 R 1 S 


icon 
l.o-z.U 


o n o o 
z.u-z.z 


Z.Z-Z.4 


p T {MeV/c) 
















75 


22.3=b 2.8 


19.0± 2.8 












125 


23.4± 2.8 


18.9± 2.0 


16.4± 2.1 


12. 5± 1.9 


17.7± 2.9 






175 




15. 9± 2.0 


13.8± 1.6 


13.5± 1.9 


12.4± 1.7 


17.6± 2.0 




225 




18.1± 2.9 


14.8± 1.7 


13.0± 1.6 


14.3± 2.1 


18. 1± 1.9 


25.4± 2.5 


275 




22. 0± 7.5 


12. 7± 1.7 


12. 8± 1.7 


13. 2± 1.4 


17.3± 1.6 


26.8± 2.3 


325 






10.3± 1.6 


14.0± 1.8 


12.9± 1.5 


16.8± 1.5 


24.3± 2.1 


375 






12.4± 3.2 


12. 2± 1.7 


12. 2± 1.6 


17.1± 1.6 


23.2± 2.0 


425 








10.6± 1.4 


13.4± 1.5 


15. 2± 1.5 


22. 7± 1.8 


475 








12. 3± 2.1 


13.6± 1.4 


17.0± 1.6 


20.8± 1.8 


525 








9.2± 2.4 


13. 5± 1.4 


15. 2± 1.5 


19.1± 2.2 


575 










11.8± 1.3 


14.7± 1.9 


17.6± 1.8 


625 










16.6± 2.2 


14.4± 1.5 


19.6± 1.9 


675 










14.0± 2.6 


15. 1± 1.5 


16.2± 1.9 


725 












13.0± 1.4 


14.4± 2.4 


775 












12.0± 1.3 


12. 8± 1.9 


825 












11. 7± 1.5 


12. 7± 1.4 


875 












9.7± 1.5 


12. 1± 1.9 


925 














10.4± 2.0 


975 














7.9± 1.4 


1025 














9.0± 1.5 


1075 














7.9± 1.6 


1125 














7.5± 1.1 



TABLE X. Invariant yields for 3 He in 10% most central Au+Pb collisions in units of 10 3 c 2 / GeV 2 . Errors listed are 
systematic and statistical combined in quadrature. 



Rapidity 


1.0-1.2 


1.2-1.4 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


2.4-2.6 


p T {MeV/c) 


















150 


8.6± 1.6 


7.9± 1.0 


6.7± 1.0 


8.6± 1.7 


8.0± 2.9 








250 


9.9± 4.4 


8.1± 1.0 


7.4± 0.7 


8.4± 0.8 


10.3± 1.1 


9.8± 1.3 


15.7± 3.8 




350 






8.4± 1.0 


7.4± 0.7 


10.1± 0.8 


12. 1± 1.1 


15.8± 1.5 


22. 2± 3.2 


450 






5.6± 1.6 


8.3± 0.9 


9.7± 0.8 


11. 7± 1.0 


16.8± 1.3 


19.1± 1.8 


550 








6.6± 0.8 


8.9± 0.9 


11.7± 0.9 


15.6± 1.2 


21.6± 1.8 


650 








7.6± 1.4 


8.6± 0.9 


12.4± 1.0 


14.3± 1.1 


21. 9± 1.8 


750 










8.0± 0.9 


11. 5± 1.0 


14.0± 1.1 


19.2± 1.6 


850 










7.2± 0.9 


11.9± 1.2 


14.7± 1.3 


17.2± 1.6 


950 










6.6± 1.4 


9.2± 0.9 


11. 8± 1.1 


15. 7± 1.7 


1050 












9.1± 0.9 


12. 7± 1.1 


13.9± 1.7 


1150 












10. 2± 1.1 


10.9± 1.0 


12. 9± 1.6 


1250 












8.1± 1.3 


9.9± 0.9 


11. 2± 1.3 


1350 












6.3± 1.2 


10.8± 1.0 


9.6± 1.1 


1450 














9.3± 0.9 


9.8± 1.0 


1550 














6.5± 0.8 


7.9± 0.9 


1650 














5.4± 0.9 


6.1± 0.8 


1750 














5.0± 0.9 


5.7± 0.7 


1850 
















4.4± 0.6 


1950 
















3.2± 0.5 


2050 
















2.6± 0.4 



3G 



TABLE XI. Invariant yields for 3 He in 10% to 38% most central Au+Pb collisions in units of lCT 3 c 2 /GeV 2 . Errors listed 
are systematic and statistical combined in quadrature. 



Rapidity 


l.U-l.Z 


10 1/1 
l.z-i.4 


1 A 1 « 
1.4-1.0 


1 « 1 Q 
l.u-l.o 


i q o n 
l.o-z.U 


o n o o 
z.u-z.z 


O O O A 


p T (MeV/c) 
















100 


8.2± 3.0 


6.4± 2.3 


6.0=b 1.8 


6.2± 2.6 








ouu 




ft 0_|_ 1 Q 
O.ZZt 1.(5 


c q_i_ i n 
O.oit l.U 


4.0 It U. / 


o.oit u.y 


Q Q_|_ 1 Q 


i i q_|_ o n 

ii.yit z.u 


500 








4.2± 0.9 


3.9± 0.6 


8.2± 0.9 


13.2± 1.2 


700 








4.9± 1.8 


5.0± 1.0 


7.8± 1.0 


11. 1± 1.0 


900 










3.7± 1.0 


8.8± 1.1 


7.7± 1.0 


1100 












5.4± 0.9 


7.8± 0.9 


1300 












2.9± 0.9 


6.7± 0.8 


1500 














3.9± 0.8 



TABLE XII. Invariant yields for 3 He in 38% to 66% most central Au+Pb collisions in units of W~ 3 c 2 /GeV 2 . Errors listed 
are systematic and statistical combined in quadrature. 



Rapidity 


1.0-1.2 


1.2-1.4 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


p T (MeV/c) 
















100 


4.8± 2.3 


1.9± 0.9 


2.4± 1.0 










300 






2.5± 0.6 


1.9± 0.4 


2.3± 0.5 


2.9± 0.6 


7.0± 1.4 


500 








1.5± 0.5 


2.5± 0.5 


2.9± 0.4 


5.2± 0.6 


700 










2.5± 0.6 


2.2± 0.4 


4.4± 0.5 


900 










1.4± 0.6 


2.7± 0.5 


3.6± 0.6 


1100 












1.7± 0.5 


3.5± 0.6 


1300 














2.4± 0.5 


1500 














1.7± 0.4 



TABLE XIII. Invariant yields for tritons in 10% most central Au+Pb collisions in units of 10 3 c 2 /GeV 2 . Errors listed are 
systematic and statistical combined in quadrature. 



Rapidity 


1.0-1.2 


1.2-1.4 


1.4-1.6 


1.6-1.8 


1.8-2.0 


2.0-2.2 


2.2-2.4 


Pt 
















50 


12.5+ 3.6 














150 


16.0+ 1.8 


11.8+ 1.5 


9.8+ 1.5 










250 


13.6+ 1.9 


12.2+ 1.5 


9.9+ 2.1 


11.8+ 2.9 


9.3+ 1.3 


15.9+ 2.5 




350 




13.6+ 1.7 


9.6+ 1.5 


12.2+ 2.4 


10.4+ 3.2 


16.0+ 1.8 


16.1+ 2.4 


450 




9.2+ 2.4 


10.2+ 1.5 


10.1+ 1.7 


11.7+ 2.5 


16.0+ 1.7 


18.5+ 2.1 


550 






10.6+ 1.4 


10.0+ 1.5 


10.5+ 3.1 


15.8+ 1.6 


18.7+ 1.9 


650 








9.0+ 1.1 


11.4+ 2.0 


15.3+ 1.7 


17.5+ 2.0 


750 








10.5+ 1.5 


9.3+ 1.6 


14.8+ 1.9 


17.0+ 1.9 


850 










10.7+ 1.5 


16.7+ 2.2 


14.8+ 2.0 


950 










11.2+ 1.3 


14.1+ 2.0 


13.9+ 2.5 


1050 












12.9+ 1.4 


15.7+ 1.8 


1150 














12.2+ 2.0 


1250 














14.3+ 1.9 


1350 














9.1+ 1.2 


1450 














9.7+ 1.3 
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TABLE XIV. Invariant yields for 4 He in 10% most central Au+Pt collisions in units of 10 5 c 2 /GeV 2 from data taken in 
the 1998 run of E864. Errors listed are systematic and statistical combined in quadrature. 



Rapidity 


1.4-1.0 


1 R 1 S 
l.O-l.o 


i q o n 


o n o o 
z.u-z.z 


O O O A 
Z.Z-Z.4 


p T {MeV/c) 












350 


17.3± 1.9 


17.4± 2.0 


22. 1± 2.6 






450 


18.1± 2.4 


16.7± 1.7 


23.3± 2.5 


26.0± 3.1 




550 


17.4± 2.4 


20.3± 2.6 


19.4± 1.9 


28.0± 3.3 


43.7± 11.4 


650 


13.4± 2.1 


18.0± 2.2 


21. 7± 2.3 


30.5± 3.5 


44.5± 11.3 


750 


20.4± 5.4 


17.2± 2.3 


21. 1± 2.4 


29.8± 3.5 


45. 6± 11.4 


850 




16.9± 2.3 


21. 2± 2.3 


31. 7± 3.8 


47.1± 11.8 


950 




17.9± 2.8 


20.3± 2.2 


28.9± 3.5 


45. 8± 11.4 


1050 




15. 5± 3.0 


23.6± 2.8 


27.6± 3.5 


42. 3± 10.6 


1150 






19.8± 2.3 


27.1± 3.4 


38. 5± 9.7 


1250 






18.7± 2.3 


27.6± 3.4 


41. 6± 10.6 


1350 






23.4± 4.0 


23.9± 3.0 


37.0± 9.5 


1450 








25. 8± 3.3 


32.0± 8.4 


1550 








21. 1± 2.8 


30. 7± 8.1 


1650 








22. 3± 3.3 


29. 1± 7.6 


1750 








20.8± 3.5 


21. 5± 5.8 


1850 








15.3± 2.7 


21. 2± 5.7 


1950 










18.5± 5.1 


2050 










13.4± 3.9 


2150 










16.5± 5.2 



TABLE XV. Invariant yields for 6 He and 6 Li in Au+Pt collisions in units of 10 8 c 2 /GeV 2 . Errors listed are systematic 
and statistical combined in quadrature. 



(y,p T (MeV/c)) 


6 He 


6 Li 


(1.6-1.8,500-1000) 


6.0 ± 1.8 




(1.8-2.0,0-500) 


10.5 ± 4.6 




(1.8-2.0,500-1000) 


5.2 ± 1.2 


32 ± 15 


(1.8-2.0,1000-1500) 


3.3 ± 1.2 


24.9 ± 6.1 


(2.0-2.2,500-1000) 


13.8 ± 1.2 


33.7 ± 6.9 


(2.0-2.2,1000-1500) 


10.4 ± 1.2 


37.5 ± 7.8 


(2.0-2.2,1500-2000) 




34 ± 21 



TABLE XVI. Invariant yields for 7 Li and 7 Be in Au+Pt collisions in units of 10 8 c 2 /GeV 2 . Errors listed are systematic 
and statistical combined in quadrature. 



(y,p T (GeV/c)) 


7 Li 


7 Be 


(1.6-2.2,0-1.8) 


.92 ± .4 




(2.0-2.2,0-1.8) 




1-3 ±:i 
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